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An ab initio downfolding method is formulated to construct low-dimensional models for 
correlated electrons. In addition to the band downfolding by constrained random phase ap- 
proximation formulated for 3D models, screening away from the target layer (chain) is further 
involved. Eliminating the off-target degrees of freedom, namely, dimensional downfolding yields 
ah initio low-dimensional models. The method is applied to derive a 2D model for a layered 
superconductor LaFeAsO, where the interlayer screening crucially makes the effective interac- 
tion short ranged and reduces the onsite Coulomb interactions by 10-20 % from the 3D model 
for the 5 iron-3d orbitals. 
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Discoveries of outstanding materials with low- 
dimensional (low-D) anisotropy such as cuprate^^ and 
iron-based^' superconductors have stimulated studies on 
low-D electronic models. '^-' ID or 2D simplified models 
have contributed in clarifying strong-correlation and fluc- 
tuation effects characteristic of the low dimensionality. 
However, to a large extent, the studies rely on ad hoc 
adjustable parameters as in the Hubbard models. 

To overcome the difficulty of density functional the- 
ory (DFT) in strongly correlated materials, a three-stage 
method combining ah initio DFT with high-accuracy 
solvers for lattice models has been proposed,*"^''' stim- 
ulated by fruitful physics and potential applications re- 
vealed in correlated electron systems. In this three-stage 
scheme, a global band structure calculated by DFT such 
as the local density approximation (LDA) is downfolded 
into a small number of bands called "target bands" 
near the Fermi level, for instance, by the constrained 
random phase approximation (cRPA).^^^' Electronic de- 
grees of freedom far from the Fermi level are eliminated 
in this step by the partial trace summation, leaving a 
lattice model for maximally localized Wannier orbitals 
(MLWO)^^-' for the target bands. This scheme has en- 
abled quantitative studies beyond the ad hoc parameter 
tuning, when the target-band models are solved by high- 
accuracy solvers. This ah initio downfolding has been 
successfully applied to various materials. ^"""'^"'-' 

However, so far, the downfolding derives models in 3D 
space. It does not offer a real connection to the exten- 
sively studied low-D models mentioned above, thus leav- 
ing a "missing link" . 

The purpose of this letter is to bridge the derived 3D 
ah initio models and low-D physics, thereby connect this 
missing link to open a new low-D studies from first prin- 
ciples. Our scheme is regarded as the dimensional down- 
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folding in real space, supplementing the original band 
downfolding in energy space. The scheme is general and 
works particularly well for quasi-low-D systems. 

Advantages of the reduction to low-D models lie princi- 
pally in numerical feasibility. The derived interaction for 
3D models has a long-range tail oc 1/r, because the con- 
ventional downfolding excludes metallic screenings com- 
pletely. So, even for quasi-low-D systems, we have to con- 
sider a 3D model whose layers (or chains) interact with 
each other by long-range interactions. The present di- 
mensional downfolding largely simplifies the problem not 
only to tractable low-D models but also to short-range- 
interaction models when metallic screening from layers 
(chains) other than the target one is considered. 

As an example, we apply it to derive an effective 2D 
model for LaFeAsO. We show that interlayer screenings 
reduce onsite Coulomb interactions by 10-20 % and fur- 
ther remove the long-range part of the screened interac- 
tion. Our formalism justifies a multi-band 2D Hubbard 
model for LaFeAsO from first principles. 

Below we take an example of quasi-2D systems. Our 
goal is to derive an effective low-energy 2D model as 



(T RR' n7n 
2 Z^ Z^ Z^l ^"'R-'iR'^nR'^mR'^mR' 
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+ JmR«R' «kamR'«nR«mR' +«,^R«„ROmR'«mR') j ' (1) 

where a^^^^ ("^nR) i^ ^ creation (annihilation) oper- 
ator of a target-band electron with spin a at the 
nth MLWO in the unit cell at R in the 2D layer. 
The single-particle levels and transfers are given 
by tmRriR'- Screened Coulomb and exchange inter- 
actions are t/mR„R' = (0mR0„R'|M^|0TOR0„R') and 
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jj^jO), where W is an effective electron 



Coulomb interaction for the 2D model. 

We employ cRPA in a two-step procedure to derive 
W for the 2D model. First, we follow the conventional 
band downfolding^' to derive the 3D model for the target 
bands: We divide the total polarization into two parts; 
a polarization x* within the target bands and the rest 
X^ ■ The division is well-defined when the target bands 
are isolated from other bands. In the usual cRPA for 
deriving a 3D model (3D-cRPA), x* is excluded and the 
screened interaction is calculated from Wr = v/{l — x'"'*^) 
with the bare Coulomb interaction v. 




Fig. 1. (Color online) Schematic diagram for effective interac- 
tions between electrons at ri and r2 in the target layer, screened 
by intra- and inter-layer polarizations x*(r,r'). Dotted lines are 
a bare Coulomb interaction, while solid lines with arrows are 
electron propagators. See the text for details. 



In the next cRPA step to derive a 2D model for the 
target layer (2D-cRPA), we eliminate degrees of freedom 
in the other layers after including the screening by the 
electrons in the other layers. Figure 1 (a) shows an RPA 
diagram (symbolically for the second-order one) of the ef- 
fective interaction between electrons at ri and r2 in the 
target layer while r and r' denote polarization points of 
X*. The electrons at ri and r (or r' and Y2) interact via v 
(dotted line). The screening of the intralayer interaction 
within the target layer is classified into several processes; 
the panel (b) shows the screening by a polarization in the 
other layers, while the panel (c) describes that by an in- 
terlayer polarization between the target and the other 
layers. The panel (d) is the screening by a polarization 
within the target layer itself. This process should be ex- 
cluded in the present 2D-cRPA downfolding, while (b) 
and (c) are included. 

In the first step, the 3D-cRPA becomes accurate and 
the vertex correction can be ignored, when the eliminated 
band energy in the denominator of the propagator is far 
from the Fermi level. In the second step, the 2D-cRPA 
becomes accurate if the numerators are small in the ver- 
tex corrections derived from small interlayer components 
of renormalized transfers and/or interactions. The small 
parameter in the vertex correction is roughly the ratio 



of the screened interlayer to onsite interactions obtained 
by full RPA. We do not consider it here. 

Practical calculations proceed as follows: The target- 
band polarization x* is given by 
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where {^/jak}, {-Eak}, and {/ak} arc the Bloch states, 
their energies, and occupancies, respectively, and the 
prime summation runs over the target bands only. The 
idea for the dimensional downfolding is to use the Fourier 
transform along the interlayer c axis [z grids) represent- 
ing the axis complementary to the target layer, -'^^^ 



XG.fi'Az^z'Y- 



V- eKq^-^^^)^x^'(q)e"'^'*^^'-^^^'(3) 



q±GsG' 



with q = q|| + q± and G = G|| + Gj^. Here, || (±) rep- 
resents an orientation parallel (complementary) to the 
layer. Suppose the interlayer hopping is switched off. 
Then Xq g'q (■^'•^O ^^ block diagonal in the z-z' plane. 
In quasi-2D materials, this block diagonal structure is 
still essentially retained as we see in Fig. 2, which en- 
ables excluding the polarization within the target layer 



as 



with 



XG„G;,q„ (2.^') = XG„G;,q,(2:^')Acut(^,^') (4) 



{0 (z,z' e target layer), 
1 (otherwise). 

The target layer is defined by a unit layer complementary 
to the c axis, bounded by the middle of the Bravais- 
lattice points. Then, we go back to the reciprocal space 
with the inverse Fourier transform 

X^„G;,q„(5^,.9l)=E^"^'"^^G„Gi,qJ-,-')^^^-^'. (6) 



Since the original lattice periodicity along the c axis is 
broken due to the cutting, in this transform, reciprocal 
lattice vector g^ for the superlattice is used, instead of 
q^+Gj^. The symmetric dielectric function^^' is 

XG||G;,q||(5-L,5l)+XG||G;q||(5-L,5l) 
~^,G' (qil) = kG " " ", , ^|, ; "'" (7) 

|qii +y\n\\ +y'\ 

with t/ = Gj|+gj^ and x^ being the 3D-cRPA polarization. 
Matrix elements of the screened Coulomb interaction are 



qiiSG' 



iq||R 



Pmq,i(^)2gg'(q||)ft*q„(^'), (8) 



where 1/ is a crystal volume and 
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The self-energy due to ofF-target propagators may also 
be calculated in 2D-cRPA, which modifies the disper- 
sion. However, in the same spirit with neglecting the 
vertex correction, we do not consider it. On the other 
hand, to avoid the double counting, the interaction al- 
ready considered in LDA has to be subtracted, where 
it is dominantly for the Hartree term and can easily be 
implemented. ^^^ 

Now we apply this formalism to derive a 2D model 
for LaFeAsO. Ab initio density- functional^^^ calculations 
were performed with plane-wave-basis-set code of Tokyo 
Ab initio Program Package^^' with LDA^^^ and norm- 
conserving pseudopotentials.^^-* The details for an iron 
pseudopotential are found in ref. 13. The experimental 
structure of LaFeAsO was taken from ref. 23. The cutoff 
energies in wavefunctions and charge densities were set 
to 100 Ry and 900 Ry, respectively, and a 5x5x5 fc-point 
sampling was employed. The polarization was expanded 
in plane waves with a 20-Ry cutoff and the total num- 
ber of bands considered in the polarization was set to 
130. The Brillouin-zone (BZ) integral on wavevector was 
evaluated by the generalized tetrahedron method.^'*-' 

We remark some details. Sufficient cutoff energy is nec- 
essary for correct decay of Xqq, (q) at large G — G'. Oth- 
erwise a noise in its Fourier transform vt, ^, iz.z') of 

eq. (3) brings about an incorrect depolarization in the in- 
verse transform Xq q, {g±,g'j_), obtained after the cut- 
ting [Eq. (4)] . We confirmed that a cutoff as large as 20 
Ry is sufficient in the present case. We also need a careful 
treatment of poles at £'ak-i-q = ^^k in eq. (2), for which, 
to reduce an error by discretized k sum, we utilize 

fak+q - //3k df{E) 



Ealti+q — EjSk 



dE 



S E' - Ev 



E=E' 



(10) 



with i?'=(i?Qk+q+i?^k)/2, based on the central-difference 
approximation to df{E)/dE ai E = E' . 

Figure 2 (a) shows a plot of XGiiGnqn (^> ^') in eq. (3) 
with G||=G||=q||=0. We see five bright spots, where po- 
larizations are large. The matrix is nearly block diagonal, 
indicating that the target-band polarization is confined 
in each layer. Thus, we can successfully exclude the po- 
larization formed within the target layer. The cutting 
region defined in eq. (5) is that inside the dashed line. 

We plot in Fig. 2 (b) diagonal elements of 2D-cRPA 
macroscopic dielectric function em = l/e~^ [(green) dots] 
obtained via inverse of eq. (7). We also show a full-RPA 
result [(red) crosses] obtained without cutting and a 3D- 
cRPA result [(blue) triangles] obtained by totally exclud- 
ing X* from the screening. We see diverging em in the 
q||+^-^0 limit in the full-RPA and 2D-cRPA results, in- 



dicating the metallic screening in both of the cases. 

Figure 2 (c) displays matrix elements of a 2D-cRPA 
screened Coulomb interaction, C/,„o„r in eq. (8), de- 
noted by (green) dots, for LaFeAsO as a function of 
the distance between the centers of the MLWOs; r = 
|('/>«R|r-|(/>„R) - {4>mo\r\(t>mo)\- In this plot, we exemplify 
the case m=n=dxy Intersite interactions were found to 
be orbital independent. Full-RPA [(red) crosses] and 3D- 
cRPA [(blue) triangles] results are also plotted. 

The interlayer screening generates a qualitatively dif- 
ferent feature in the model: It turns the long-ranged in- 
teraction in 3D-cRPA into short ranged. We find that for 
LaFeAsO the interactions larger than 0.1 eV are limited 
in the region up to the next-nearest neighbors (r < 4A), 
implying the screening length ~ AA determined by the 
interlayer channel. It justifies the use of a 2D short- 
ranged-interaction model from first principles. 

A careful analysis is needed to reach the convergence 
for the number of screening layers n. Taking five kj_ 
points in BZ, for instance, simulates the number of 
screening layers n = 2 each above and below the tar- 
get layer in a supercell. Since we need the interaction 
parameter in the limit n — ^ 00, the data is extrapolated 
by U{n) = U{oo) - [[/(oo) - U{0)]e-^"/'^^ . Here, a repre- 
sents the effective screening length in the layer unit and 
t/(0) is the 3D-CRPA value. The inset of Fig. 2 (c) shows 
the fitting of the onsite intraorbital Uxy We obtain U{oo) 
= 2.75 eV and a = 2.1. The interlayer screening reduces 
the onsite U by 10 to 20 % from the 3D-cRPA values. 

Table I summarizes the onsite interaction parameters 
derived from 3D-cRPA and extrapolated 2D-cRPA. The 
diagonal interactions are compared in the upper two 
rows. The interlayer screening is nearly orbital indepen- 
dent and U2D is uniformly reduced by '^0.4 eV from U^d. 
This is because the screening length is much larger than 
the spread of MLWOs. In contrast to the diagonal part 
of the interaction, the interlayer screening does not affect 
off-diagonal exchange interactions (lower two rows). 

Recent variational Monte Carlo simulation^*^ has re- 
vealed that the use of the 3D-cRPA parameters'^-* results 
in an overestimate of the antifcrromagnetic ordered mo- 
ment as ^-^2.0 jiB in contrast to the experimental value 
0.3-0.6 /iB-^'^'^^'' A part of the discrepancy is solved by 
considering La-4/-elcctron screening ignored in the 3D- 
cRPA calculation in ref. 11, yielding 10% reduction of 
U}^^ However, the present 2D-cRPA is crucial in fur- 
ther reducing 10-20 % of the U values, thus leading to 
a nearly perfect agreement with the experiment. An ap- 
plication to K-ET conductors also indicating that this 
method offers a promising versatile scheme will be dis- 
cussed elsewhere. The present method opens up a way 
to study low-dimensional physics from first principles. 

In summary, we have formulated a general ab initio 
scheme to derive an effective low-D model by a down- 
folding of a 3D model. The formalism eliminates the de- 
grees of freedom for layers (chains) other than the target 
one, leaving an ab initio low-D model. After decompos- 
ing the polarization in the target band into layer/chain- 
by- layer /chain contributions in the real space, the low-D 
model is obtained with the interaction screened by the 
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Fig. 2. (Color online) (a) Xq g q (■^'■^') (^^ a.u.) in the z-z' plane with G|i=G||=q|i=0. Inside the dashed line is the cutting region 
defined in eq. (5). (b) Dielectric functions as a function of Iqii +G\ obtained by fuU-RPA [(red) crosses], 2D-cRPA [(green) dots], and 
3D-cRPA [(blue) triangles], (c) Screened Coulomb interactions of LaFeAsO as a function of distance r between centers of MLWOs. 
Only the interactions between d^y MLWOs, Uxy{i')=Uj;yOxyfi, are plotted. The intersite interactions arc nearly the same between 
the inter- and intra-orbital pairs. The same symbols as (b) are used for interaction values. The 3D-cRPA interaction is scaled as 
"'^/[^M*'^^ (^~*''^)''] (^"li"^ curve). (Inset) Convergence of onsite intraorbital Uxy{r = 0) as a function of number of screening layers n 
to n — >■ oo. Solid curve is the fitted U{n) (see the text). 



Table I. Comparison between onsite 3D-cRPA and 2D-cRPA interaction parameters for eq. (1). Upper (lower) two rows give diagonal 



(off-diagonal) interaction. The unit is eV. In the column labels, (1,2,3,4,5) := {xy,yz,3z 
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intcrlaycr/chain channel from RPA polarizations away 
from the target layer/chain. For metallic systems such as 
LaFeAsO, this screening is essential because it deletes the 
long-range part of the interactions. These results justify 
the 2D short-ranged models as effective ab initio models. 
This scheme enables first-principles studies of strongly 
correlated materials with low-D anisotropies heretofore 
extensively studied only by models with ad hoc parame- 
ters. 
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